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FINITE  DEFORMATION  EFFECTS  IN  PLASTICITY  ANALYSIS* 


E*  H*  I^ft 

Department  of  Mechanical  Engineering 
Aeronautical  Engineering  and  Mechanics 
Rensselaer  Polytechnic  Institute,  Troy,  New  York  12181 


ABST^^ 

The  non-linear  kinematics  of  the  cordiination  of  elastic  and 
plastic  deformations  at  finite  strain  provides  the  mathematical 
structure  to  examine  aspects  of  elastic-plastic  analysis  more 
succinctly  than  is  possible  with  the  approach  based  on  infinitesimal 
elastic  strain. 

Kinematic  hardening  represents  the  anisotropic  component  of 
strain  hardening  by  a  back  stress  ct.  Application  of  current  theory  for 
finite  deformation  Incorporates  the  effect  of  finite  rotation  by 
using  the  Jaumann  derivative  in  the  evolution  equation  for  a. 

This  approach  predicts  oscillating  shear  stress  for  monotonlcally 
increasing  simple  shear  strain  but  this  anomaly  can  be  eliminated 
by  adopting  a  physically  more  meaningful  modified  Jaumann  derivative. 

1.  INTRODUCTION 

Structural  metals  can  often  be  deformed  to  large  strains  with¬ 
out  fracturing.  The  rational  design  of  many  engineering  processes 
and  structures  demands  stress  and  deformation  analysis  in  the  presence 
of  finite  strain.  In  order  to  anticipate  and  hence  prevent  the  gen¬ 
eration  of  cracks,  high  residual  stresses  or  other  forming  defects  in 
the  manufacture  of  structural  components,  stress  analysis  must  apply 
throughout  a  body  and  so  elastic-plastic  theory  must  be  utilized. 

The  neglect  of  the  small  elastic  strain  compared  with  plastic  strain, 
and  hence  the  adoption  of  rigid-plastic  theory,  prevents  the  determ¬ 
ination  of  stresses  in  the  rigid  regions  which  may  well  comprise  most 
of  the  work-piece  in  a  metal  forming  process  for  example. 


*For  presentation  at  the  International  Conference  on  Numerical  Methods 
in  Industrial  Forming  Processes  at  University  College  of  Swansea, 

U.K. ,  July  1982. 


The  early  development  of  elastic-plastic  theory  was  based  on 
infinitesimal  deformation  theory  on  the  basis  of  which  the  total 
strain  was  equal  to  the  sum  of  elastic  and  plastic  components 
with  a  similar  summation  law  applying  also  for  strain  rates.  At 
finite  strains  there  is  a  coupling  between  elastic  and  plastic 
deformation  since  plastic  flow  occurs  in  a  material  already 
stressed  to  yield,  and  hence  subject  to  elastic  strains,  and 
these  two  components  interact  in  the  nonlinear  kinematics  of 
finite-deformation  theory.  In  the  currently  commonly  used 
approach  to  finite-deformation  analysis,  the  summation  of  elastic 
and  plastic  strain  rates  to  give  the  total  strain  rate  is  adop¬ 
ted.  The  significance  of  this  assumption  is  examined  in  the 
light  of  the  nonlinear  klnematical  theory.  The  latter  gives  a 
precision  to  the  kinematics  which  permits  many  aspects  of  the 
theory  to  be  Investigated  more  succinctly. 

2.  ELASTIC-PLASTIC  KINEMATICS 

A  body  is  considered  in  its  initial  undeformed  unstressed 
configuration  labelled  X  in  Fig.  1.  After  a  deformation  which 


Fig.  1.  Elastic-plastic  deformation. 


involves  both  elastic  and  plastic  strain,  the  body  occupies  the 
configuration  labelled  x,  the  body  having  been  deformed  at  time  t 
according  to  the  motion 

X  *■  x(X,t)  (2.1) 

In  order  to  uncouple  the  elastic  and  plastic  deformations,  the 
body  is  considered  to  be  destressed  from  the  configuration  jc  to 
the  configuration  p.  It  is  assumed  that  the  destressing  to  zero 
stress  is  purely  elastic,  although  the  theory  can  be  modified  to 
include  materials  which  exhibit  a  strong  Bauschinger  effect 
involving  plastic  flow  during  the  destresslng  to  zero  stress  [1]. 
Since  the  stress  is  zero  in  the  configuration  g,  the  elastic 
strain  is  zero  and  hence  p  exhibits  the  plastic  strain.  Tills  is 
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also  the  plastic  strain  present  In  the  configuration  x  since  only 
the  elastic  strain  has  changed  during  the  destress IngT 


The  elastic-plastic  deformation  which  has  taken  place  during 
the  motion  (2.1)  Is  expressed  by  the  deformation  gradient  F 


F  -  Bx/BX  (F^j  -  3x^/3Xj) 


(2.2) 


The  mapping  X  p  with  deformation  gradient  F^  expresses  the 
plastic  defonsation  In  both  the  configuration  P  and  x.  The 
mapping  p  x  with  the  deformation  gradient  F^  expresses  the 
elastic  ^formation  In  the  configuration  x.  Because  removing  the 
surface  tractions  from  a  body  which  has  been  subjected  to  non- 
homogeneous  plastic  strain  leaves  It  In  a  state  of  residual 
stress,  destresslng  Involves  considering  the  body  sectioned  Into 


vanishingly  small  elements  and  thus  the  mapping  to  configuration 
p  from  either  X  or  x  Is  not  differentiable.  F^  and  F^  are  then 
^Int  functions  of  position  but  not  deformation  gradients.  This 
causes  no  difficulty  since  the  objective  Is  to  develop  a  con¬ 
stitutive  relation  for  the  stress  in  the  configuration  x  and  the 
corresponding  mapping  from  X  to  x  is  differentiable  so  that  the 
deformation  gradient  F  does'exlst.  Considering  the  body  sec¬ 
tioned  Into  small  elements  to  achieve  destresslng  to  zero  stress 
Is  analogous  to  machining  away  parts  of  a  specimen  to  measure 
residual  stresses.  Such  partition  does  not  affect  the  basic 
elastic  law  so  that  using  this  In  combination  with  measurements 
of  change  In  strain,  the  stresses  in  the  original  whole  component 
can  be  deduced.  In  practice,  test  specimens  subjected  to  homo¬ 
geneous  strain  are  commonly  used  and  then  unloading  the  specimen 
does  Indeed  destress  It. 


The  configurations  X,  x  and  p  can  be  obtained  and  measured 
experimentally,  and  the  mpplng  x'  -»■  x  Is  Identical  geometrically 
with  the  resultant  of  the  sequence  X  ^  p  x.  The  chain  rule 
Chen  generates  the  relation 


F  -  fV  (2.3) 

The  generally  non-commutatlve  structure  of  this  matrix  product 
relation  Indicates  that  It  Is  not  compatible  with  the  summation 
law  for  elastic  and  plastic  strain  or  strain  rate.  It  should 
perhaps  be  pointed  out  that  the  sequence  of  maps  X-^(X,C)-»x(X,C) 
for  varying  time  t  cannot  be  carried  out  physically,  because^  for 
example,  plastic  flow  cannot  be  generated  at  zero  stress,  but 
since  the  configurations  can  be  achieved  and  then  measured  by 
passing  through  x,  this  does  not  affect  the  applicability  of  the 
relation  (2.3). 

The  unstressed  state  p  Is  not  unique  since  arbitrary  rota¬ 
tion  leaves  it  unstressed.'  Thus,  without  loss  of  generality, 
unstresslng  and  hence  the  Inverse,  elastic  deformation  f  ,  can 
be  chosen  to  be  pure  deformation  without  rotation  and  thus  can 


b«  expressed  by  the  synnetric  defonnatloa  grsdient  matrix  V* 
This  choice  simplifies  some  parts  of  the  analysis. 
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Substituting  (2.3)  Into  the  expression  for  total  finite 
Lagrange  strain 


E  -  (F^F  -  I)/2 


(2.4) 


(superscript  T  denoting  transpose  and  I  the  unit  matrix)  ex¬ 
presses  the  total  strain  In  terms  of  the  elastic  and  plastic 
strains  calculated  from  F^  and  F^  by  relations  analogous  to 
(2.4)! 


fP  e*  F**  + 


(2.5) 


For  large  plastic  strains  this  Is  clearly  far  removed  from 
summation  of  elastic  and  plastic  strains. 

In  plasticity  theory  strain  rate  Is  appropriately  defined 
In  terms  of  the  velocity  field  (y(x,t)  In  the  current  elastlcally- 
plastlcally  deformed  configuration's)  as  the  symmetric  part  of 
the  velocity  gradient  L: 

L  -  3v/3x  -  Ov/3X)(ax/3x)  -  F  F"^  (2.6) 

where  F  expresses  a  time  derivative  at  fixed  X.  that  Is  at  a 
fixed  raterlal  particle.  Expressing  L  as  the  sum  of  Its  sym¬ 
metric  and  antl-symmetrlc  parts 


L  -  D  +  W  (2.7) 

determines  the  deformation  rate  or  stretching  tensor  D  and  the 
spin  W.  The  former  expresses  the  rate  of  strain  about  the  cur¬ 
rent  configuration  x  which  is  appropriate  for  plasticity  an¬ 
alysis.  Plasticity' theory  Is  commonly  termed  Incremental  or 
flow  type,  more  akin  to  fluid  behavior  than  to  solid,  the  initial 
configuration  X  playing  a  minor  role  unlike  the  case  of  elas¬ 
ticity  theory.' 

Substituting  (2.3)  with  V®  into  (2.6),  incorporating  (2.7), 
yields 

p  -  p®  +  y®pPy®  ig  +  y®wPy®  |g  (2.8) 

where  p®,  W®,  p^  and  are  the  symmetric  and  antl-symmetrlc 

parts  of  V®V®  and  F^F^  respectively  and  subscripts  S  and  A 
denote  the  symmetric  and  antl-symmetrlc  parts  respectively. 

Now  y  Is  the  symmetric  part  of  (I  -f  3u  /3g)  where  u  "x-p 
Is  the  elastic  displacement.  Since  elastic  strains  are'us^lly 
of  the  order  (yield  stress  divided  by  elastic  modulus)  *-10~^, 
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V*  is  close  to  Che  unit  nacrlx  I.  There  are  exceptions  to  this 
circumstance,  for  example  in  explosive  generated  shock  waves  in 
metals  when  volumetric  elastic  strain  can  be  of  the  order  unity. 
When  'v  I,  (2.8)  approximates  the  strain  rate  summation  re- 
elation 

p  -  D*  +  D**  (2.9) 

since  is  synaetric  and  anti-sysnetric .  Thus  the  common 
assumption  in  finite-deformation  elastic-plastic  theory  usually 
provides  a  good  approximation  to  the  theory  based  on  nonlinear 
klnesiatics.  However,  in  considering  certain  aspects  of  the 
structure  of  the  theory,  the  precision  of  the  nonlinear  kln- 
ematlcal  theory  permits  a  more  incisive  investigation. 

3.  ELASTIC-PLASTIC  CONSTITUTIVE  RELATIONS 

Plasticity  exhibits  a  "flow"  or  "incremental"  type  of  re¬ 
sponse  to  stress  in  that  the  plastic  strain  rate  (or  some  char¬ 
acteristic  of  it)  is  determined  by  the  stress  variation,  so  that 
the  plastic  strain  is  determined  by  Integrating  the  rate  re¬ 
lations  through  the  stress  history.  The  final  deformation  thus 
depends  on  the  stress  history  and  similarly  the  stress  depends  on 
the  deformation  history.  For  a  strain  hardening  material  the 
strain-rate  is  given  in  terms  of  the  current  state  (stress  and 
hardening  parameters)  and  the  stress  rate.  For  deformation  of 
most  structural  metals  at  ambient  temperatures,  the  plasticity 
law  is  of  first  order  in  time  rates  on  each  side  of  the  con¬ 
stitutive  equation  so  that  any  monotonically  increasing  function 
of  time  can  replace  the  time  variable  without  modifying  the 
resulting  stress  and  strain  histories.  For  example  a  monoto¬ 
nically  Increasing  displacement  or  strain  variable  would  be 
appropriate.  Materials  which  exhibit  such  a  response  are  termed 
time-rate  independent  though  they  are  of  flow  or  incremental 
type.  This  paper  is  concerned  with  such  materials,  although  the 
finite-deformation  kinematics  discussed  applies  equally  well  to 
rate  dependent  laws  which  usually  apply  at  higher  temperatures. 

Since  the  plasticity  law  involves  plastic  strain  rate,  the 
formulation  of  elastic-plastic  theory  requires  that  the  plas¬ 
ticity  law  and  elasticity  law  be  substituted  into  the  total 
strain  rate  relation  (2.8),  or  the  approximation  to  it  (2.9),  to 
provide  an  expression  for  the  resultant  strain  rate  D  in  terms  of 
stress,  stress  rate  and  variables  which  express  the  influence  of 
the  previous  history  of  the  stress  or  deformation. 

Since  elasticity  is  usually  expressed  as  a  function  relation 
for  the  stress  in  terms  of  the  deformation,  a  rate  form  of  this 
relation  must  be  derived  for  substitution  into  (2.8)  or  (2.9). 

In  the  usual  approach  consideration  is  limited  to  small  elastic 
strains  and  Hooke's  law  is  used  in  the  form 


(3.1) 


e  1  +  V  V  . 

®±j  E  "  E  ■'kk^lj 


friiere  e  Is  che  Infinitesimal  strain  and  t  the  Klrchhoff  stress 
[2,3].'  A  rate  form  of  this  Is  commonly' taken  to  be 


1  +  V 


ij 


VO  . 

E  “^kk^lj 


(3.2) 


0  £ 
where  t  Is  a  time  derivative  of  the  stress  tensor  and  D  cor¬ 
responds  to  the  total  rate  of  deformation  tensor  defined  In 
(2.7)  when  the  elastic  velocity  gradient  Is  used. 


Because  the  final  constitutive  relation  must  be  valid  for 
large  plastic  deformation.  It  must  be  so  for  large  total  strain 
and  rotation.  It  Is  thus  Important  to  check  that  the  constit¬ 
utive  relation  Is  objective,  i.e.  that  If  a  time-dependent  rigid 
body  motion  Is  superimposed  on  the  deformed  configuration  x,  the 
effect  on  the  stress  at  time  t  must  be  simply  to  rotate  It'ac- 
cordlng  to  the  value  of  the  superimposed  rotation  at  that  time. 

If  (3.2)  Is  substituted  Into  (2.9)  and  combined  with  an  objective 
plasticity  law  It  Is  clear  that  the  stress  rate  t  must  transform 
under  rigid  body  rotation  expressed  by  the  proper-orthogonal 
matrix  Q  (t)  as  does  D 


D  9  p  q'^  (3.3) 

There  Is  an  Infinity  of  such  rates  to  choose  from  [4] ,  but  the 
Infinitesimal  approximation  built  Into  (3.2)  prevents  a  deductive 
choice.  In  current  finite-element  elastic-plastic  programs  two 
such  rates  have  been  adopted  for  flnlte-deformatlon  analysis,  the 
Jaumann  rate  and  the  Truesdell  rate. 

In  the  transformation  of  (3.1)  to  rate  form  (3.2),  (3.1)  was 
not  formally  differentiated  because  It  Is  not  a  valid  equality  In 
the  context  of  finite  deformation  theory.  Rather  a  derivative  of 
a  type  not  uniquely  selected  was  applied  to  the  right-hand  side, 
and  the  velocity  field  was  Introduced  with  which  to  express  the 
’’rate  of  strain"  for  the  left-hand  side.  In  contrast  the  non¬ 
linear  klnematlcal  relation  (2.8)  can  be  combined  with  the  flnlte- 
deformatlon-valld  rubber  elasticity  law  to  permit  formal  dif¬ 
ferentiation  of  both  sides  of  the  elastic  constitutive  relation 
to  provide  a  strain  rate  term  which  can  be  substituted  Into 
(2.8).  This  procedure  was  carried  out  In  [5]  using  the  Jaumann 
derivative.  Since  the  same  operator  Is  applied  to  both  sides 
of  the  elastic  constitutive  relation,  an  equality  will  result. 

The  trick  Is  to  choose  the  operator  so  that  the  resulting 
strain-rate  like  term  can  be  substituted  Into  (2.8).  It  Is 
necessary  that  the  plasticity  law  will  also  mesh  conveniently 
Into  (2.8).  These  components  all  meshed  together  In  [5]  to 
yield  an  elastic-plastic  constitutive  relation  and  a 


corresponding  variational  principle.  A  simplification  of  the 
final  result  valid  for  small  elastic  strain  yielded  the  con¬ 
ventional  expression. 
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The  plasticity  analysis  In  [5]  was  based  on  Isotropic 
hardening  with  Mlses  yield  condition.  The  finite  deformation 
elasticity  lav  was  also  Isotropic  and  not  affected  by  plastic 
flow.  In  these  circumstances  the  appropriate  choice  for  the 
stress  rate  term  has  already  been  considered  by  Prager  [4] . 


The  plasticity  law  takes  the  form 


p  .  ifii  o 

^ij  h  3t..  3t  ^mn 
^  Ij  mn 


(3.4) 


so  that  zero  stress  rate  will  Imply  zero  plastic  strain  rate.  And 
yet  for  the  named  objective  stress  rate  definitions  which  Prager 
considered,  other  than  the  Jaumann  rate,  vanishing  of  the  stress 
rate  did  not  Imply  stationary  behavior  of  the  stress  Invariants. 
Thus  for  rate  definitions  other  than  the  Jaumann  rate,  the  stress 
Invariants,  and  hence  the  yield  condition,  would  change  when  no 
plastic  flow  was  taking  place.  This  is  certainly  not  consistent 
with  Isotropic  plasticity  theory.  Hence  of  the  stress  rate  de¬ 
finitions  considered.  Including  the  Truesdell  rate,  only  the 
Jaumann  rate  Is  appropriate  for  elastic-plastic  theory.  Prager 
used  the  word  "preferable"  which  seems  to  be  too  weak  a  statement. 


Other  aspects  of  elastic-plastic  theory  can  be  more  pre¬ 
cisely  Investigated  using  the  nonlinear  kinematics.  For  example 
the  last  term  In  (2.8)  represents  a  rate  of  change  of  elastic 
strain  which  would  appear  as  residual  strain  after  a  stress 
Increment  has  been  applied  and  then  removed  [6].  This  fact  was 
utilized  In  [5]  by  coupling  It  to  to  obtain  a  total  elastic 
term  which  meshed  conveniently  with  the  elastic-plastic  Incre¬ 
mental  constitutive  relation  development  already  described. 

It  was  shown  In  [1]  that  this  elastic  strain  rate,  associ¬ 
ated  with  the  occurrence  of  the  spin  while  the  body  Is  under 
stress.  Is  parallel  to  the  tangent  to  the  yield  surface  In  the 
nlne-dinenslonal  stress  and  strain  space.  This  vector  Is  thus 
normal  to  the  stress  vector  and  hence  no  power  Is  expended  by 
this  strain  rate.  An  analogous  situation  has  a  bearing  on  the 
Investigation  of  the  stability  of  plastic  flow. 

Analysis  of  the  stability  of  plastic  flow  by  Drucker  [7], 
amplified  by  Palmer,  Haler  and  Drucker  18]  considers  the  power 
(rate  of  work)  expended  In  plastic  flow.  This  rate  of  work  takes 
place  In  the  configuration  x  and  Is  expressed,  according  to  (2.8) 
by 


vhere  tr  indicates  the  trace.  Since  t  is  symmetric  this  is 
equal  to 

tr(ty'D'’y*  )  (3  5, 

For  the  case  under  consideration  associated  with  the  Mises  yield 
condition,  the  matrices  t,  and  V®  all  have  the  same  principal 
directions  so  that  the  matrix  products  in  (3.6)  are  commutative, 
hence  (3.6)  reduces  to 

tr(TDP)  (3.7) 

Thus  the  plastic  rate  of  working  by  the  stress  on  the  motion  in 
configuration  x  is  the  same  as  that  for  the  purely  plastically 
deformed  configuration  p.  This  uncouples  the  Influence  of  elastic 
strain  on  the  analysis  of  this  question,  a  significant  insight 
provided  by  the  adoption  of  nonlinear  finite  deformation  kine¬ 
matics.  Note  that  the  commutativity  of  the  deformation  associ¬ 
ated  matrices  in  (3.6)  determines  the  equality 

y®p^y®  -  p**  (3.8) 

so  that  normality  of  the  plastic  strain  race  vector  applies  also 
to  its  manifestation  in  Che  x  configuration.  This  result  depends 
on  the  adoption  of  the  Mises  yield  condition.  For  an  isotropic 
yield  condition  involving  the  stress  deviacor  invariant  J3 
as  well  as  J2,  for  example^the  Tresca  condition,  is  not  par¬ 
allel  to  T,  nor  hence  to  y®  and  (3.8)  is  not  valid. 

4.  STRESS  ANALYSIS  FOR  KINEMATIC  HARDENING 

So  far  Che  discussion  has  been  limited  to  isotropic  har¬ 
dening.  However  it  is  known  Chat  the  anisotropic  Bauschinger 
effect  becomes  more  pronounced  as  the  strain  Increases  and  can 
have  a  significant  effect  particularly  if  destresslng  occurs 
followed  by  reverse  stressing.  This  can  arise  as  Che  product 
emerges  from  a  metal-forming  process.  It  is  therefore  important 
to  incorporate  anisotropic  hardening  and  the  most  utilized  ap¬ 
proach  has  been  kinematic  hardening  or  a  combination  of  this  with 
isotropic  hardening.  In  Che  Workshop  on  Plasticity  at  Finite 
Deformation  held  at  Stanford  University  in  1981,  NagCegaal  and  de 
Jong  presented  a  paper  [9]  in  vdilch  Che  stresses  generated  by 
simple  shear  Co  large  deformation  in  materials  which  exhibit 
anisotropic  hardening  were  evaluated  using  analysis  currently 
accepted  as  valid  for  finite  deformation.  They  obtained  the 
unexpected  result,  for  a  material  which  strain  hardens  mono- 
tonlcally  in  tension,  that  Che  shear  traction  grows  to  a  maximum 
value  at  a  shear  strain  y  of  order  unity  and  then  oscillates 
with  increasing  strain.  This  paper  has  stimulated  further  study 
of  Che  commonly  accepted  approach  to  the  analysis  of  finite 
deformation  in  the  presence  of  anisotropic  hardening  and  a  sug¬ 
gested  modification  of  the  theory  follows. 
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In  conformity  with  current  practice  for  finite  deformation 
analysis  In  the  presence  of  kinematic  hardening,  [9]  adopted  an 
evolution  equation  for  the  back  stress  or  shift  tensor  a  (the 
current  center  of  the  yield  surface)  which  relates  the  Jaumann 
derivative  of  a  to  the  plastic  strain  rate.  This  Incorporates 
effects  of  finite  rotation  and  ensures  objectivity  of  the  evo¬ 
lution  equation  under  superimposed  rigid-body  rotations. 

Fig.  2  illustrates  the  kinematics  of  simple  shear  at  finite 
strain.  The  velocity  field  Is  given  by 

Vj^  -  kx2.  V2  -  v^  -  0  (4.1) 

« 

k  being  the  rate  of  shear  strain  y.  The  material  spin  W  Is 
uniform  and  constant  and  represents  an  angular  velocity  of  mag¬ 
nitude  k/2.  The  evolution  equation  takes  the  form 

a  -  CD^  (4.2) 


Fig.  2.  Simple  shear  In-  the  direction. 

where  C  is  a  scalar  function.  Thus  the  material  derivative  which 
expresses  the  rate  of  change  of  o  relative  to  the  fixed  coordi¬ 
nates  shown  in  Fig.  2  is  given  by 

a  -  cpP  +  Wo  -aW  (4.3) 

It  is  pointed  out  in  [10]  that  in  view  of  the  constant  spin 
matrix  W  the  spin  terms  in  (4.3)  cause  a  to  rotate  continuously 
and  that  this  causes  the  oscillations  in  stress  already  referred 
to. 


It  Is  shown  in  [10]  and  evident  from  Fig.  2  that  no  loci  of 

•  material  particles  ever  rotate  through  more  than  a  limited  angle 

as  t  For  example  the  material  particles  lying  initially  on 

*  the  X2  axis  OA^  only  approach  rotation  through  ir/2  radians  as 

c  -)>a>.  IhLs  poses  a  paradox  since  the  back  stress  a  is  generated  by 
micro-mechanisms  embedded  in  the  material  and  it  seems  implaus¬ 
ible  that  Of  would  rotate  continuously  if  the  physical  source  of 
its  effects  did  not.  The  mlcro-mechanlsms  could  be  dislocations 
on  crystallographic  slip  planes  piled  up  against  grain  boundaries 
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or  Inclusions  or  residual  stresses  In  crystallites  generated  by 
the  Inhooogeneous  polycrystalllne  structure. 


Since  the  back  stress  a  expresses  an  average  over  the  many 
active  mlcro-mechanlsms.  It  seems  reasonable  to  select  an  av¬ 
erage  direction  embedded  In  the  material  the  rotation  of  which 
constitutes  a  contribution  to  the  rotation  of  a.  A  further 
contribution  arises  from  the  mechanisms  being  activated.  The 
latter  Is  represented  by  the  first  term  on  the  right-hand  side  of 
(4.3).  Study  of  the  micro-mechanics  of  plastic  flow  In  an  aniso¬ 
tropic  medium  will  be  required  to  determine  the  orientation  of 
the  relevant  material  loci,  but  the  macroscopic  kinematic  har¬ 
dening  law  provides  evidence  from  which  this  could  be  Inferred. 

In  [10]  we  suggested  that  the  eigen-vector  corresponding  to  the 
eigen-value  of  a  of  maximum  absolute  value  provides  a  plausible 
choice.  The  rotation  of  th^  locus  of  material  elements  so  ori¬ 
ented  determines  the  spin  W  .  A  modified  Jaumann  derivative 
based  on  this  spin  Is  defined: 

a*a-Wa  +  aW  (4.4) 


and  the  corresponding  evolution  equation  for  a  becomes 


(4.5) 


giving 


a  -  CpP  +  W*a  -  aW*  (4.6) 

It  Is  shown  in  [10]  that  the  modified  evolution  equation  (4.4) 
Is  objective. 


In  formulating  the  plastic  flow  relations  care  must  be  taken 
since  the  yield  condition  will  be  rotating  with  the  tensor  a 
and  the  effect  of  this  must  be  eliminated  from  the  definition  of 
the  stress  rate  In  the  expression  for  the  plastic  strain  rate. 
This  Is  treated  in  a  manner  analogous  to  that  for  the  evolution 
equation  by  defining  a  modified  Jaumann  derivative  of  the  stress 
based  on  the  rotation  of  the  a  tensor.  The  summation  of  strain 
rates  (2.9)  was  adopted  slnce^  as  already  discussed,  it  provides 
a  satisfactory  approximation. 

Calculations  were  carried  out  [10]  for  kinematic  hardening 
although  the  theory  developed  applies  for  combined  kinematic- 
isotropic  hardening.  A  model  involving  linear  hardening  In 
tension  was  assumed  with  constants  chosen  to  approximate  an 
aluminum  alloy. 

Fig.  3  shows  the  deduced  stress  variation  based  on  rigid- 
plastic  analysis.  For  this  problem  of  homogeneous  deformation 
rlgld-plastlc  analysis  Is  satisfactory  since  no  rigid  regions 
occur  In  the  body.  Since  rlgld-plastlc  theory  Involves  In¬ 
compressible  deformation  (which  is  consistent  with  the  velocity 
boundary  conditions)  the  stress  can  only  be  determined  to  with- 
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Fig.  3.  Stress  versus  shear  strain  for  rlgld-plastlc  model. 

in  an  arbitrary  uniform  hydrostatic  stress.  The  normal  stress 
plotted,  Is  In  fact  the  stress-devlator .  With  boundary 

velocities  prescribed,  elastic-plastic  analysis  determines  a 
unique  stress  field  which  Is  devlatoric  and  agrees  with  the 
rlgld-plastlc  solution. 

Fig.  3  shows  that  the  method  In  current  use  based  on  the 
conventional  Jaumann  derivative  predicts  the  stress  oscillation 
in  both  shear  and  normal  components  obtained  in  [9].  Isotropic 
hardening  with  linear  hardening  In  tension  yields  linear  har¬ 
dening  in  shear  with  no  normal  stresses  generated  The  new  ap¬ 
proach  determines  smooth  hardening  curves  in  shear  and  normal 
stress.  No  oscillations  occur  since  now  the  hardening  parameter 
a  rotates  towards  a  limiting  position  as  t  increases.  The  new 
approach  gives  close  agreement  with  the  other  solutions  at  small 
strains  and  a  gradual  reduction  of  the  shear  hardening.  The 
principal  direction  of  the  anisotropic  hardening  rotates  towards 
the  xi  axis  and  this  causes  the  reduced  hardening  in  the  shear 
curve. 

Comparison  of  the  analyses  with  the  conventional  and  the 
modified  Jaumann  derivatives  indicates  that  up  to  strains  near 
0.5  the  difference  in  the  solutions  is  small.  For  strains  near 
2  the  difference  reaches  about  40%  and  grows  rapidly  with  in¬ 
creasing  strain.  These  results  suggest  that  current  codes  for 
stress  analysis  at  large  strains  which  incorporate  kinematic 
hardening  should  be  reassessed  if  the  possibility  of  serious 
errors  in  stress  evaluations  is  to  be  avoided. 
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The  non-linear  kinematics  of  the  combination  of  elastic  and  plastic  defor¬ 
mations  at  finite  strain  provides  the  mathematical  structure  to  examine  aspects  of 
elastic-olastic  analysis  more  succinctly  than  is  possible  with  the  approach  based 
on  Infinitesimal  elastic  strain. 

Kinematic  hardening  represents  the  anisotropic  component  of  strain  hardening 
by  a  back  stress  g.  Application,  of  current  theory  for  finite  deformation  incorpor¬ 
ates  the  effect  of  finite  rotation  by  using  the  Jaumann  derivative  in  the  evolution 
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equation  for  a.  This  approach  predicts  oscillating  shear  stress  for 
monotonlcally  increasing  simple  shear  strain  but  this  anomaly  can  be 
eliminated  by  adopting  a  physically  more  meaningful  modified  Jaumsnn 
derivative. 
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